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These notes give examples of how suitably defined geometrical objects encode in their fractal 
structure thermal critical behavior. The emphasis is on the two-dimensional Potts model for which 
two types of spin clusters can be defined. Whereas the Fortuin-Kasteleyn clusters describe the 
standard critical behavior, the geometrical clusters describe the tricritical behavior that arises when 
including vacant sites in the pure Potts model. Other phase transitions that allow for a geomet- 
rical description discussed in these notes include the superfluid phase transition and Bose-Einstein 
condensation. 
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I. INTRODUCTION 



The quest for understanding phase transitions in terms 
of geometrical objects has a long history. One of the ear- 
lier examples, due to Onsager, concerns the superfluid 
phase transition in liquid 4 He — the so-called A transition. 
During the discussion of a paper presented by Gorter, 
Onsager [l[ made the following remark: "As a possible in- 
terpretation of the A-point, we can understand that when 
the concentration of vortices reaches the point where they 
form a connected tangle throughout the liquid, then the 
liquid becomes normal." Feynman also worked on this 
approach and summarized the idea as follows [2j|: "The 
superfluid is pierced through and through with vortex 
line. We are describing the disorder of Helium I." This 
approach focuses on vortex loops, i.e., one-dimensional 
geometrical objects, which form a fluctuating vortex tan- 
gle. As the critical temperature T\ is approached from 
below, the vortex loops proliferate and thereby disorder 
the superfluid state, causing the system to revert to the 
normal state. The A transition is thus characterized by 
a fundamental change in the typical vortex loop size. 
Whereas in the superfluid phase only a few small loops 
are present, close to T\ loops of all sizes appear. The sud- 
den appearance of arbitrarily large geometrical objects is 
reminiscent of what happens in percolation phenomena 
at the percolation threshold where clusters proliferate. 
Even on an infinite lattice, a percolating cluster can be 
found spanning the lattice. 

A second example, due to Feynman is related to 
Bose-Einstein condensation. Here, the relevant geomet- 
rical objects are worldlines. In the imaginary-time for- 
malism, used to describe quantum systems at finite tem- 
perature T, the time dimension becomes compactified, 
t = —it, with < r < h/k^T, where fee is Boltzmann's 
constant. Because of periodic boundary conditions, the 
worldlines then form closed loops. At high temperatures, 
where the system behaves more or less classically, the 
individual particles form separate closed loops wrapping 
only once around the imaginary time axis. Upon lowering 
the temperature, these small loops, describing single par- 
ticles, hook up to form larger exchange rings. A particle 
in such a composite ring (see Fig. [1} moves in imaginary 



time along a trajectory that does not end at its own start- 
ing position, but ends at that of another particle. Hence, 
although the initial and final configurations are identical, 
the particles in a composite ring are cyclically permuted 
and thus become indistinguishable [3| . Figure [1] gives 
an example of three particles, labeled 1,2, and 3. After 
wrapping once around the imaginary time axis particle 1 
ends at the starting position of particle 2, which in turn 
ends after one turn around the imaginary time axis at 
the starting position of particle 3. That particle, finally, 
ends at the starting position of particle 1. In this way, 
the three particles are cyclically permuted, forming the 
cycle (1,2,3). Being part of a single loop which winds 
three times around the imaginary time axis, the particles 
cannot be distinguished any longer. At the critical tem- 
perature, worldlines proliferate and — again as in percola- 
tion phenomena — loops wrapping arbitrary many times 
around the imaginary time axis appear, signaling the on- 
set of Bose-Einstein condensation 0, Q ■ This approach 
has been turned into a powerful Monte Carlo method by 
Ceperley and Pollock Q that can even handle strongly 
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FIG. 1: The worldlines of three particles that, after moving a 
time t = h/ksT in the imaginary time direction, are cyclically 
permuted (left panel) . The three separate worldlines can also 
be represented by a single worldline that winds three times 
around the imaginary time axis (right panel). (After Ref. |J|.) 
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FIG. 2: Snapshots of typical spin configurations of the Ising 
model on a square lattice of linear size L — 100 in the normal, 
hot phase at /3 = 0.5 /3 C (left panel) and just above the Curie 
point at /3 = 0.98 /3c (right panel). A spin up is denoted by a 
black square, while a spin down is denoted by a white one. 

interacting systems like superfluid 4 He (see Ref. Q for a 
review). 

A third example concerns the phase transition in sim- 
ple magnets. The most elementary model describing such 
a transition is provided by the Ising model, obtained by 
assigning a spin that can point either up or down to each 
lattice site. Figure [2] shows typical spin configurations for 
a square lattice in the normal, hot phase and just above 
the Curie point. For convenience, a spin up is denoted 
by a black square, while a spin down is denoted by a 
white one. From these snapshots, the relevant geomet- 
rical objects appear to be clusters of nearest neighbor 
spins in the same spin state (in the following, we will 
qualify this statement). The normal, disordered phase 
consists of many small clusters. As the Curie point T c is 
approached from above, larger clusters appear, which at 
T c start to proliferate — as in percolation phenomena. In 
the absence of an applied magnetic field, the percolating 
cluster can consist of either up or down spins, both having 
equal probability to form the majority spin state. Since 
the percolating spin clusters have a fractal structure, it is 
tempting to ask whether this structure encodes the stan- 
dard thermodynamic critical behavior, as in percolation 
theory? More generally, we wish to address in these notes 
the question: Can suitably defined geometrical objects 
encode in their fractal structure the standard critical be- 
havior of the system under consideration? To highlight 
the basic features, we consider simple models, such as the 
Ising, the Potts, and the XY model. Moreover, we study 
them mostly in two dimensions (2D) since many analyti- 
cal predictions, obtained by using Coulomb gas methods 
and conformal field theory, are available there. 

The rest of these notes is organized as follows. In the 
next section, the 2D critical Potts model is discussed. 
Central to the discussion is the equivalent geometrical 
representation of this spin model in terms of so-called 
Fortuin-Kasteleyn clusters @. The fractal structure of 
these stochastic clusters and the way the thermal criti- 



cal behavior of the Potts model can be extracted from 
it are studied in detail. In Sec. IIIIl the tricritical Potts 
model is discussed. The clusters encoding the tricriti- 
cal behavior turn out to be the naive clusters of nearest 
neighbor spins in the same spin state, which feature in 
Fig. [2 Their fractal structure is connected via a dual 
map to that of the Fortuin-Kasteleyn clusters, which en- 
code the thermodynamic critical behavior. In Sec. IIV1 
the boundaries of both cluster types are studied. The 
notes end with a summary of the main results and an 
outlook to other applications. 

II. CRITICAL POTTS MODEL 

A. Fortuin-Kasteleyn Representation 

The Potts model is one of the well studied spin mod- 
els in statistical physics @. It is defined by consider- 
ing a lattice with each lattice site given a spin variable 
Si = 1, 2, • • • , Q that can take Q different values. In its 
standard form, the spins interact only with their nearest 
neighbors specified by the Hamiltonian 

H = -K y £(8 at ,. j -1), (1) 

m 

where K denotes the coupling constant. Nearest neigh- 
bor spins notice each other only when both are in the 
same spin state, as indicated by the Kronecker delta. 
The Potts model is of particular interest to us as for 
Q = 2 it is equivalent to the Ising model, while in the 
limit Q — ► 1 it describes ordinary, uncorrelated percola- 
tion. The notation Y]uj\ is to indicate that the double 
sum over the lattice sites, labeled by i and j, extends over 
nearest neighbors only. The partition function Z can be 
written as 

Z = Tr e-^ = Tr J] [(1 - p) + p 5 Si<Sj ] , (2) 
(«> 

where denotes the inverse temperature, and the trace 
Tr stands for the sum over all possible spin configura- 
tions. In writing Eq. ([2]), use is made of the identity 

e^°i-^ = (l-p)+p6 Si , Sj , (3) 

with p = 1 — c~P , where here and in the sequel we set 
the coupling constant K to unity. The identity §5§ can be 
pictured as setting bonds with probability p/[(l—p)+p] = 
p between two nearest neighbor spins in the same spin 
state for which 5 SuSj = 1. When two nearest neighbor 
spins are not in the same spin state, S Si ^ Sj — 0, then with 
probability (1 — — p) = 1 the bond is not set, i.e., 

never. It thus follows, that the partition function can be 
equivalently written as 

z FK = J2p b (i-p)~ b+a Q Nc , (4) 

{r} 
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where {r} denotes the set of bond configurations realized 
on a total of B bonds on the lattice. A given configu- 
ration is specified by b set and b not set bonds between 
nearest neighbor spins in the same spin state, and a pairs 
of nearest neighbor spins not in the same spin state (for 
which the bonds are never set). Together they add up to 
the total number of bonds, B = b + b + a, so that the 
exponent b + a in Eq. ^ can also be written as B — b. 
Only spins connected by set bonds form a cluster. The 
exponent Nq in Eq. ^ denotes the number of clusters, 
including isolated sites, contained in the bond configura- 
tion F. The factor Q Nc arises because a given cluster can 
be in any of the Q possible spin states. Equation ((4]) is 
the celebrated Fortuin-Kasteleyn (FK) representation of 
the Potts model @ . It gives an equivalent representation 
of that spin model in terms of FK clusters obtained from 
the naive geometrical clusters of nearest neighbor spins 
in the same spin state, discussed in the Introduction, by 
putting bonds with a probability p — 1 — between 
nearest neighbors. As geometrical clusters are split up 
in the process, the resulting FK clusters are generally 
smaller and more loosely connected than the geometrical 
ones. 

Not only does the FK representation provide a geo- 
metrical description of the phase transition in the Potts 
model, it also forms the basis of efficient Monte Carlo al- 
gorithms by Swendsen and Wang [loj , and by Wolff [ll[ , 
in which not individual spins are updated, but entire FK 
clusters. The main advantage of the nonlocal cluster up- 
date over a local spin update, like Metropolis or heat 
bath, is that it substantially reduces the critical slowing 
down near the critical point. 



B. FK Clusters 

The results of standard percolation theory [l2j also ap- 
ply to FK clusters. In particular, the distribution l n of 
FK clusters, giving the average number density of clus- 
ters of mass n, takes near the critical point the asymp- 
totic form 

e n ~n- T e~ en . (5) 



we start by considering the radius of gyration R n , 

i=l i,j=l 

with Xj the position vectors of the sites and x = 
(1/n) Y^h=i x i the center of mass of the cluster. Asymp- 
totically, the average (R n ) scales with the cluster mass n 
as 

(R n ) ~ n 1 / 15 , (8) 

which defines the Hausdorff, or fractal dimension D. The 
average radius of gyration (Rn) gives the typical linear 
size of a cluster of mass n. A second length scale is pro- 
vided by the correlation length £, which diverges close 
to T c with an exponent v as £ ~ \T — T c \~ v . Both are 
related via 

(Rn)=ZR(n6), (9) 

where R is a scaling function, cf. Eq. ([5|). From the 
asymptotic behavior ((Sj> , the divergence of the correla- 
tion length, and the vanishing © of the parameter 9 as 
T c is approached, the relation 



follows, connecting the critical exponent v to the fractal 
dimension D of the clusters and a. 

The fractal dimension can also be related to the en- 
tropy exponent r as follows. At criticality, the mass n 
of a cluster is distributed over a volume of typical linear 
size (Rn), so that 

nl n ~ l/(Rn) d , (11) 

with d the dimension of the lattice. This leads to the 
well-known expression 

r=^ + l, (12) 

in terms of which the correlation length exponent reads 
v=(r-l)/da. 



The first factor, characterized by the exponent r, is an 
entropy factor, measuring the number of ways a cluster of 
mass n can be embedded in the lattice. The second factor 
is a Boltzmann weight which suppresses large clusters 
when the parameter 9 is finite. Clusters proliferate and 
percolate the lattice when 9 tends to zero. The vanishing 
is characterized by a second exponent a defined via 

9 oc \T-T C | 1/a . (6) 

As in percolation theory [k|, the values of the two ex- 
ponents specifying the cluster distribution uniquely de- 
termine the critical exponents. To obtain these relations, 



C. Improved Estimators 

To see how physical observables, such as the magneti- 
zation m and the magnetic susceptibility x are repre- 
sented in terms of FK clusters, we consider the Ising 
model in the standard notation with the spin variable 
Si = ±1 for simplicity. 

The correlation function (S^Sj) has a particular sim- 
ple representation. When the two spins belong to two 
different FK clusters 

<$Sj> = ^ E ^ = (13) 
Si,Sj=±l 
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while when they belong to the same cluster 

(SiSj) = l £ S^ = l. (14) 

That is, if Ci denotes the FK cluster to which the spin 
Si belongs and Cj the one to which Sj belongs, then 

(S i S j ) = S Ci ,c r (15) 
For the susceptibility \ = J2ij(SiSj) in the normal 



phase, Eq. (fTS]) gives 



ij 



{C} 



(16) 



where the sum J2{c} * s over a ^ FK dusters, and ric 
denotes the mass of a given cluster. In terms of the FK 
cluster distribution £ n , the susceptibility can be written 
as 



(17) 



Note that in percolation theory E2|, the ratio 
J2 n n2 ^n/ J2n n ^n denotes the average cluster size. Since 
in the Ising model all L d spins are part of some FK clus- 
ter, we have the constraint 



1. 



(18) 



It thus follows that the right hand of Eq. (fTT)) precisely 
gives the average size of FK clusters. In other words, 
this geometrical observable directly measures the mag- 
netic susceptibility of the Ising model. From the asymp- 
totic form ([5]), and the divergence x ~ \T — T c |~ 7 of the 
susceptibility when the critical point is approached, the 
relation 7 = (3 — t)/<t between the critical exponent 7 
and the cluster exponents a and r follows. 

Also the magnetization m has a simple geometrical rep- 
resentation [1 21 ] - In an applied magnetic field H, a spin 
cluster of mass nc has a probability oc exp(/3nc H) to 
be oriented along the field direction, and a probability 
cx exp(—(3ncH) to be oriented against the field direc- 
tion. The difference between these probabilities gives the 
magnetization mc per spin in the cluster, 



toc = tanh(/3?ic-ff). 



(19) 



Close to the critical temperature and in the thermody- 
namic limit L d — > co, the largest cluster dominates, and 
tanh(/?n max _ff) — > ±1 for this cluster, depending on its 
orientation. The magnetization of the entire system (per 
spin) then becomes 

m = ±P OQ , (20) 

where = n max / L d gives the fraction of spins in the 
largest cluster — the so-called percolation strength. Be- 
cause of the constraint (p~8|) . it is related to the FK cluster 
distribution via 



1 



£ 



nir, 



(21) 



where the prime on the sum indicates that the largest FK 
spin cluster is to be excluded. The magnetization van- 
ishes near the critical point asm~ \T-T c \ p . Together 
with the asymptotic behavior of the cluster distribution, 
Eq. ([20]) with Eq. flU) gives the relation = (r - 2) /a. 

These geometrical observables (average cluster size and 
percolation strength) are called improved estimators be- 
cause they usually have a smaller standard deviation than 
the spin observables. 

The results just derived for the Ising model also apply 
to the rest of the critical Potts models [8[ . In this way, the 
thermal critical exponents of these models are completely 
determined by the exponents a and r, characterizing the 
FK cluster distribution. Specifically, 



a = 2 - 



rj = 2 + d 



t-1 



a 

t — 3 

7^1' 



T-2 

a 
t-1 

da ' 



3 - r 



7 = 



(22) 



as in percolation theory [121 ] . The exponent rj, determin- 
ing the algebraic decay of the correlation function at the 
critical point, is related to the fractal dimension via 



Consequently 



D = 2 — 77). 



j/u = 2D-d. 



D. Critical Exponents 



(23) 



(24) 



The critical exponents of the 2D Q-state Potts model 
are known exactly [l3| . It is convenient to parametrize 
the models as 



= —2 cos(7t/k), 



(25) 



with 2 > R > 1. For the Ising model (Q = 2) k = 4/3, 
while for uncorrelated percolation (Q — ► 1) R = 3/2. 
The correlation length exponent v and the exponent 77 
are given in this representation by fl3l ]: 



1 



2/T,l 



V 



(26) 



where 2/t,i is the leading thermal exponent. The next- 
to-leading thermal exponent ?/t,2 reads j/t,2 = 4(1 — k), 
which is negative for R > 1, implying that the corre- 
sponding operator is an irrelevant perturbation. The 
other critical exponents can be obtained through stan- 
dard scaling relations. The parameter R is related to the 
central charge c, defining the universality class, via [l4| 



c = 1 



6(1 -R) 2 



(27) 



Finally, the fractal dimension D of FK clusters is given 
by El El 



1 3 
D = 1+ 2R- + 8 K ' 



(28) 



5 




FIG. 3: Distribution £ n normalized to the volume L 2 of FK 
clusters in the 2D Ising model at criticality on a square lattice 
of linear size L — 512. Statistical error bars are omitted 
from the data points for clarity. The straight line is a one- 
parameter fit through the data points with (minus) the slope 
fixed to the predicted value r = 31/15 = 2.06667. The fit 
illustrates that asymptotically the distribution is algebraic, 
as expected at criticality. 



FIG. 4: Log-log plot of the percolation strength Poo of geo- 
metrical and FK clusters at criticality in the 2D Ising model 
as a function of the linear system size L. The straight lines 
0.988281 L' - 0527 for geometrical and 1.00558 L" 1248 for FK 
clusters are obtained from two-parameter fits through the 
data points. Statistical error bars are smaller than the symbol 
sizes. 



which gives D = 15/8 for the Ising model and D = 91/48 
for uncorrelated percolation. 

To demonstrate that FK clusters actually percolate at 
the critical point, Fig.[3]shows the distribution t n of these 
clusters in the 2D Ising model at criticality (0 = 0) on 
a square lattice of linear size L = 512. With D = 15/8, 
it follows from Eq. (fT2"]l that the entropy exponent takes 
the value r = 31/15. The straight line, obtained through 
a one-parameter fit with the slope fixed to the predicted 
value, shows that asymptotically the FK cluster distri- 
bution has the expected behavior. 

E. Swendsen-Wang Cluster Update 

The theoretical predictions (|2"2")1 can be directly verified 
through Monte Carlo simulations, using the Swendsen- 
Wang cluster update [l0|. Instead of single spins, entire 
FK clusters are considered units to be flipped as a whole 
in this approach. Standard finite-size scaling theory ap- 
plied to the percolation strength P M and the average clus- 
ter size x gives the scaling laws 

P 00 =L-P/"P(L/Q, x = L?/"X(L/0, (29) 

with P and X scaling functions. Precisely at T c , these 
scaling relations imply an algebraic dependence on the 
system size L, allowing for a determination of the expo- 
nent ratios f3/v (see Fig. 0} and j/v. Using these geo- 
metrical observables as improved estimators for the mag- 
netization and susceptibility, respectively, we arrived at 
the estimates for the Ising model (Q = 2) [l7| 

/3/u = 0.1248(8) « 1/8, 

7/1/ = 1.7505(12) « 7/4, (30) 

where the right hands give the known values for the Ising 
critical exponents. These estimates illustrate first of all 



that FK clusters indeed encode the thermal critical be- 
havior of the Ising model. Moreover, they also illustrate 
that measuring geometrical observables gives excellent 
results for the critical exponents. The data were fitted 
over the range L = 64 — 512, using the least-squares 
Marquardt-Levenberg algorithm. 

In Ref . [IH , the fractal dimension of FK clusters were 
obtained from analyzing their distribution. This method 
gives less accurate results than applying finite-size scaling 
to improved estimators. The main problem is related to 
the fitting window. The fitting range cannot be started 
at too small cluster sizes, where the distribution has not 
taken on its asymptotic form yet, while too large clus- 
ter sizes, which are generated only a few times during a 
complete Monte Carlo run, are also to be excluded be- 
cause of the noise in the data and finite-size effects. The 
results depend sensitively on the precise choice of the fit- 
ting window. 

F. Geometrical Clusters 

Figure [3] suggests that the geometrical spin clusters 
also percolate right at the Curie point of the Ising model. 
To demonstrate this to be the case, Fig. [5] shows the dis- 
tribution of these clusters at criticality. Asymptotically, 
the distribution indeed shows algebraic behavior, imply- 
ing that clusters of all size appear in the system. It is 
therefore natural to investigate the exponents associated 
with the percolation strength P G (see Fig. H]) and the av- 
erage size x G °f these geometrical clusters. (The super- 
script "G" refers to geometrical clusters.) Using finite-size 
scaling, as for the FK clusters, we arrived at the estimates 

m 

j3 G jv = 0.0527(4) w 5/96 = 0.0521 • • ■ , 

j G /v = 1.8951(5) w 91/48 = 1.8958 •■ • . (31) 
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FIG. 5: Distribution normalized to the volume L 2 of ge- 
ometrical clusters in the 2D Ising model at criticality on a 
square lattice of linear size L — 512. Statistical error bars 
are omitted from the data points for clarity. The straight line 
is a one-parameter fit through the data points with (minus) 
the slope fixed to the value r G = 379/187 = 2.02674. The fit 
illustrates that asymptotically the distribution is algebraic at 
criticality. 



In obtaining these estimates we included percolating clus- 
ters. When excluding them, as was done in Ref. [l9|], the 
estimates become less accurate [17J]. The entropy expo- 
nent which follows from these results is r G = 379/187, 
corresponding to the fractal dimension D G — 187/96. 

It should be stressed that only in 2D geometrical clus- 
ters percolate right at the critical temperature. In higher 
dimensions, geometrical clusters percolate in general too 
early at a lower temperature, and their fractal structure 
is unrelated to any thermodynamic singularity. 

The 2D exponents ([31)) are not related to the criti- 
cal behavior of the Ising model and the question arises: 
What do these exponents describe? 



III. TRICRITICAL POTTS MODEL 

A. Dual Map 

When the pure 2D Potts model is extended to include 
vacant sites, it displays in addition to critical also tricriti- 
cal behavior at the same critical temperature T c [20( . The 
tricritical behavior is known to be intimately connected 
to the critical behavior, and both critical points share 
the same central charge. To demonstrate this connec- 
tion, note that for a given c, Eq. (|2"7|) yields two solutions 
for R: 



K± 



13 - c ± y/{c-2b){c- 1) 



12 



(32) 



with = 1, where R = R+ > 1 and hence R— < 1. 

Stated alternatively, the substitution R with l/R leaves 
the central charge (|2T|) unchanged, c(R) = c(l/R). When 
applied to the parametrization ([23)1 of the critical Potts 
branch, this so-called dual map yields the parametriza- 



/ <3 t 



-2 cos(7tk) 



(33) 



of the tricritical branch (the superscript "t" refers to the 
tricritical point). Various results for the critical point 
[lU can be simply transcribed to the tricritical point by 
using this dual map, leading to [13, [U 

1 . „ 3 ( „ _ 3 



V 



(34) 



2k' '' 4k' 
while the next-to-leading thermal exponent becomes 

Z/T,2=4-i (35) 

To preserve relation (|23[) under the dual map, the fractal 
dimensions of the geometrical and FK clusters must also 
be related by the map R — » l/R [H|, HH . This gives 

+ (36) 

which is indeed the correct fractal dimension of geomet- 
rical clusters [H], HH|. In other words, the geometrical 
clusters can, as far as their scaling behavior is concerned, 
be considered shadows of the FK clusters. The use of the 
word "shadow" will become clear when we consider the 
cluster boundaries in the next section. 



B. Ising & its Q l = 1 Potts Shadow 

Equation ([33)) gives as fractal dimension of the geomet- 
rical clusters of the Ising model [R = 4/3) D G — 187/96, 
implying via Eq. ([T2")l r = 379/187, in accordance with 
what we found numerically [T7[. Note that with R = 4/3, 
Eq. ([33)) gives Q t = 1. That is, the tricritical model de- 
scribed by the geometrical clusters of the Ising model is 
the diluted Q t = 1 Potts model. Both models share the 
same central charge c — 1/2. 

The alert reader may have noticed a curiosity con- 
cerning the thermal exponents. According to Eq. ([3~4")l . 
the correlation length exponent v t takes the value v l = 
I/j/t.i = 8 / 15 in thc diluted Q* = 1 Potts model 
(R = 4/3). Yet, in our numerical investigation [ItJ of 
the geometrical clusters of the Ising model, we seem to 
observe the correlation length exponent v = 1 of the 
Ising model. Hence, v and not the tricritical exponent 
v l appears in Eq. |[3"T)) . In fact, what we see is the tri- 



critical next-to-leading thermal exponent (j35p . which for 
the diluted Q t — 1 Potts model happens to take the same 
value as the leading thermal exponent of the Ising model, 
Ut,i = 1 for R = 4/3. 



Vt,2 



IV. HULLS & EXTERNAL PERIMETERS 
A. FK Clusters 

When clusters percolate at a certain threshold, their 
boundaries necessarily do too. In the context of uncorre- 
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Fractal Dimensions 



(a) 



FIG. 6: In both panels, a piece of the same single FK cluster 
of nearest neighbor sites (filled circles) connected by bonds 
(black links) is shown. Two different external boundaries can 
be defined: (a) The hull (dark filled circles) is found by allow- 
ing a random walker tracing out the boundary to move only 
over set bonds, (b) The external perimeter (dark filled circles) 
is found by allowing the random walker to move to a near- 
est neighbor on the cluster boundary irrespective of whether 
the connecting bond is set or not. The external perimeter, 
which contains two sites less than the hull for this boundary 
segment, is therefore a smoother version of the hull. 



lated percolation in 2D, external cluster boundaries can 
be traced out by a biased random walker as follows [2(| . 
The algorithm starts by identifying two endpoints on a 
given cluster, and putting the random walker at the lower 
cndpoint. The walker is instructed to first attempt to 
move to its nearest neighbor to the left. If that site is 
vacant, the walker should try to move straight ahead. If 
that site is also vacant, the walker should try to move to 
its right. Finally, if also that site is vacant, the walker 
is instructed to return to the previous site, to discard 
the direction already explored, and to investigate the (at 
most two) remaining directions in the same order. When 
turning left or right, the walker changes its orientation 
accordingly. The procedure is repeated iteratively un- 
til the upper endpoint is reached. The other half of the 
boundary is obtained by repeating the entire algorithm 
for a random walker instructed to first attempt to move 
to its right rather than to its left. 



For FK clusters, being built from bonds between near- 
est neighbor sites with their spin in the same spin state, 
one can imagine two different external boundaries (see 
Fig. O- First, one can allow the random walker to move 
along the FK boundary only via set bonds. This defines 
the hull of the cluster. Second, one can allow the ran- 
dom walker to move to a nearest neighbor site on the FK 
boundary irrespective of whether the bond is set or not. 
This defines the external perimeter of the cluster, which 
is a smoother version of the hull. Since boundaries are 
clusters themselves, they too are characterized by a frac- 
tal dimension and a distribution like Eq. ([5]) specified by 
two exponents a and r. 



The fractal dimensions of the hulls (H) and external 
perimeters (EP) of FK clusters are given by [H, H3, HH 



1 
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EP 



1 



1 

2R' 



(37) 



As for clusters, the average hull and external perimeter 
sizes diverge at the percolation threshold. Let 7h and 
7ep denote the corresponding exponents, then because 
of Eq. (HU) with d = 2 and Eq. §7$) 



Jr/v = h, jep/v=1/k, 



(38) 



where a single correlation length exponent v is assumed. 

For illustrating purposes, Fig.[7]shows the distribution 
of the two boundaries of FK clusters in the Ising model 
at criticality. The straight lines are one-parameter fits 
through the data points with the slopes fixed to the ex- 
pected values. Although the estimates for Dr and Dep 
we obtained, using finite-size scaling applied to the im- 
proved estimators at criticality, are compatible with the 
theoretical conjectures [22j |, the achieved precision is less 
than the one we reached for the clusters themselves. The 
reason for this is as follows. While including percolat- 
ing clusters when considering the mass of the clusters, 
we ignore them in tracing out cluster boundaries. Be- 
cause of the finite lattice size, large percolating clusters 
have anomalous small (external) boundaries, so that in- 
cluding them would distort the boundary distributions. 
Moreover, the Grossman- Aharony algorithm [26| used to 
trace out cluster boundaries generally fails on a percolat- 
ing cluster as its boundary not necessarily forms a single 
closed loop any longer. However, as we explicitly demon- 
strated for the cluster mass 17J , disregarding percolating 
clusters leads to strong corrections to scaling, and there- 
fore to less accurate results. 



C. Geometrical Clusters 

For geometrical clusters, where the bond between near- 
est neighbor sites with their spin in the same spin state 
is so to speak always set, hulls and external perimeters 
cannot be distinguished, and 



D 



EP- 



(39) 



The fractal dimension of the boundary is gotten from 
that of the hull (|3T[) of FK clusters by applying the dual 
map k — ► l/^i yielding 29] 



(40) 



Since FK clusters have two boundaries, while geometri- 
cal clusters have only one, geometrical clusters have less 
structure and can be considered shadows of FK clusters 
under the dual map, as far as their scaling behavior is 
concerned. 
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FK clusters 




FIG. 7: Distribution normalized to the volume L 2 of the hulls 
and external perimeters of FK clusters in the 2D Ising model 
at criticality on a square lattice of linear size L = 512. Statis- 
tical error bars are omitted from the data points for clarity. 
The straight lines are one-parameter fits through the data 
points with the slopes fixed to the expected values. For clar- 
ity, the external perimeters distribution is shifted downward 
by two decades. 

Geometical clusters 




FIG. 8: Distribution normalized to the volume L 2 of the hulls 
geometrical clusters in the 2D Ising model at criticality on a 
square lattice of linear size L = 512. Statistical error bars 
are omitted from the data points for clarity. The straight line 
is a one-parameter fit through the data points with the slope 
fixed to the expected value. 



Again for illustrating purposes, Fig. [8] shows the dis- 
tribution of the hulls of geometrical clusters in the Ising 
model at criticality. The slow approach to the asymptotic 
form, with the associated strong corrections to scaling we 
observed for the hulls of geometrical clusters, stands out 
clearly from the other distributions. The reason for this 
is that geometrical clusters have a larger extent than FK 
clusters. On a finite lattice, percolating clusters gulp 
up smaller ones reached by crossing lattice boundaries. 
For geometrical clusters this happens more often than 
for FK clusters, so that disregarding percolating clusters 
when tracing out cluster boundaries has a more profound 
effect. In particular, the average hull size is underesti- 
mated. With increasing lattice size, the effect becomes 
smaller, as we checked explicitly [17J. 

This behavior is different from what we found in an- 
other numerical study of the hulls of geometrical clusters 



in the 2D Ising model [22|. In that Monte Carlo study, 
we used a plaquette update to directly simulate the hulls. 
Although the largest hull was omitted in each measure- 
ment, we found corrections to scaling to be virtually ab- 
sent (see Fig. 11 of that paper). This allowed us to obtain 
a precise estimate for the fractal dimension on relatively 
small lattices. 



V. CONCLUSIONS & OUTLOOK 

As illustrated in these notes, for the 2D Potts mod- 
els it is well established that suitably defined geometrical 
objects encode in their fractal structure critical behav- 
ior. In fact, two types of spin clusters exist, viz., FK and 
geometrical clusters, which both proliferate precisely at 
the thermal critical point. As emphasized before, this is 
special to 2D. In general, geometrical clusters percolate 
at an inverse temperature /3 P > [3 C . The fractal structure 
of FK clusters encodes the critical exponents of the criti- 
cal Potts model, while that of geometrical clusters in 2D 
encodes those of the tricritical Potts model. The fractal 
structure of the two cluster types as well as the two fixed 
points are closely related, being connected by the dual 
map k — > 1/R. This map conserves the central charge, 
so that both fixed points share the same central charge. 
The geometrical clusters can, as far as scaling properties 
are concerned, be considered shadows of the FK clusters. 

Up to now we considered external boundaries of spin 
clusters as clusters themselves, which necessarily perco- 
late when the spin clusters do. An alternative way of 
looking at these boundaries is to consider them as loops. 
In this approach, it is natural to extend the Ising model 
in another way and to consider the 0(A) spin models, 
with —2 < N < 2. The high-temperature (HT) represen- 
tation of the critical 0(A) spin model naturally defines a 
loop gas, corresponding to a diagrammatic expansion of 
the partition function in terms of closed graphs along the 
bonds on the underlying lattice [3(3] • The loops percolate 
right at the critical temperature, and similar arguments 
as given in these notes for spin clusters show that the 
fractal structure of these geometrical objects encode im- 
portant information concerning the thermal critical 0(A) 
behavior [3l|, 1321. This connection was first established 
by de Gennes [33| for self-avoiding walks, which are de- 
scribed by the 0(A) model in the limit A — > 0. One 
aspect in which lines differ from spin clusters is that they 
can be open or closed. It is well known from the work 
on self-avoiding walks that the loop distribution itself is 
not sufficient to establish the critical behavior, as has re- 
cently also been emphasized in Ref. [3ij . For this, also 
the total number z n = J^. z n (x.i,Xj) of open graphs of 
n steps starting at and ending at an arbitrary site Xj 
is needed. Its asymptotic behavior close to the critical 
temperature, cf. Eq. (O with Eq. ([T2"j) . 



(41) 



provides an additional exponent which together with 
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the loop distribution exponents is needed to specify the 
full set of critical exponents [32|. In Eq. (f4"Tj) , D denotes 
the fractal dimension of the closed graphs. Note that for 
spin clusters, the notion of open or closed does not apply, 
so that the analog of the exponent i9 is absent there. 

Remarkably, the HT graphs of a given critical O(N) 
model represent at the same time the hulls of the ge- 
ometrical clusters in the Q-state Potts model with the 
same central charge [22|, [2|| [H[ . To close the circle, we 
note that, as in the Potts model, including vacancies in 
the 0(N) model gives rise to also tricritical behavior. 
The tricritical point corresponds to the point where the 
HT graphs collapse. In the context of self-avoiding walks 
(N — > 0), this point is known as the point. Using the 
duality discussed in Sec. IIII Ai we recently conjectured 
that the tricritical HT graphs at the same time represent 
the hulls of the FK clusters of the Potts model with the 
same central charge c as the tricritical O(N) model [3l| . 
This connection allowed us to predict the magnetic scal- 
ing dimension of the 0(N) tricritical model, in excellent 
agreement with recent high-precision Monte Carlo data 
in the range < c < 0.7 [3a |. 

We started these notes mentioning the A transition in 
liquid 4 He in terms of vortex proliferation. In closing, we 
wish to give the present status of that picture as estab- 
lished in a very recent high-precision Monte Carlo study 
of the 3D complex |</>| 4 theory describing the transition 
[37l ]. An important observable is the total vortex line 
density v. By means of standard finite-size scaling anal- 
ysis of the associated susceptibility \ = L 3 ({v 2 ) — (v) 2 ), 
the inverse critical temperature /3 C was estimated and 
shown to be perfectly consistent with the estimate of 
a previous study directly in terms of the original vari- 
ables [38]. Unfortunately, when considering percolation 
observables, such as whether a vortex loop percolates the 
lattice, slight but statistically significant deviations from 
/3 C were found. For all observables considered, the per- 
colation threshold j3 p > (3 C . That is, from these observ- 
ables one would conclude that the vortices proliferate too 
early at a temperature below the critical one (as do ge- 
ometrical clusters in 3D). Yet, when taking the perco- 
lation threshold as an adjustable parameter, reasonable 
estimates were obtained from percolation observables for 
the critical exponents v and /3, consistent with those of 
the XY model. The problem with the percolation thresh- 
old is quite possibly related to the way vortex loops are 
traced out. When two vortex segments enter a unit cell, 
it is not clear how to connect them with the two outgo- 



ing segments. A popular choice is to randomly connect 
them, but it might well be that the resulting network is 
too extended and consequently percolates too early. It 
is in our mind conceivable that a proper prescription for 
connecting vortex segments could lead to a vortex perco- 
lation threshold right at the critical temperature, in the 
spirit of the FK construction. 

As a final remark, we note that even in cases where 
no thermodynamic phase transition takes place, the no- 
tion of vortex proliferation can be useful in understand- 
ing the phase structure of the system under considera- 
tion. An example is provided by the 3D Abelian Higgs 
lattice model with compact gauge field [3^|. In addition 
to vortices, the compact model also features magnetic 
monopoles as topological defects. It is well established 
that in the London limit, where the amplitude of the 
Higgs field is frozen, it is always possible to move from 
the Higgs region into the confined region without en- 
countering thermodynamic singularities [40| . Neverthe- 
less, the susceptibility data for various observables define 
a precisely located phase boundary. Namely, for suffi- 
ciently large lattices, the maxima of the susceptibilities 
at the phase boundary do not show any finite-size scaling. 
Moreover, the susceptibility data obtained on different 
lattice sizes collapse onto single curves without rescaling, 
indicating that the infinite-volume limit is reached. In 
Ref. [39[ it was argued that this phase boundary marks 
the location where the vortices proliferate. A well-defined 
and precisely located phase boundary across which geo- 
metrical objects proliferate, yet thermodynamic quanti- 
ties remain nonsingular has become known as a Kertesz 
line. Such a line was first introduced in the context of 
the Ising model in the presence of an applied magnetic 
field [if . 
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